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[57] ABSTRACT 

A GPS attitude receiver for determining the attitude of a 
moving vehicle in conjunction with a first, a second, a third, 
and a fourth antenna mounted to the moving vehicle. Each 
of the antennas receives a plurality of GPS signals that each 
include a carrier component. For each of the carrier com- 
ponents of the received GPS signals there is an integer 
ambiguity associated with the first and fourth antennas, an 
integer ambiguity associated with second and fourth anten- 
nas, and an integer ambiguity associated with the third and 
fourth antennas. The GPS attitude receiver measures phase 
values for the carrier components of the GPS signals 
received from each of the antennas at a plurality of mea- 
surement epochs during an initialization period and at a 
measurement epoch after the initialization period. In 
response to the phase values measured at the measurement 
epochs during the initialization period, the GPS attitude 
receiver computes integer ambiguity resolution values rep- 
resenting resolution of the integer ambiguities. Then, in 
response to the computed integer ambiguity resolution val- 
ues and the phase value measured at the measurement epoch 
after the initialization period, it computes values defining the 
attitude of the moving vehicle at the measurement epoch 
after the initialization period. 

28 Claims, 8 Drawing Sheets 
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SYSTEM AND METHOD FOR GENERATING 
ATTITUDE DETERMINATIONS USING GPS 

The United States Government may have certain rights in 
any patent issuing from this application pursuant to the 5 
Grant NAS8-36125 between NASA and Stanford Univer- 
sity. 

This application is a continuation-in-part of U.S. patent 
application Ser. No. 08/036,319 filed Mar. 24, 1993, now 
abandoned. 10 


FIELD OF THE INVENTION 

The present invention relates generally to systems and 15 
methods for generating precise attitude determinations for 
moving vehicles. In particular, it pertains to a system and 
method for generating attitude determinations for moving 
vehicles, including aircraft and spacecraft, using the Global 
Positioning System (GPS). 20 

BACKGROUND OF THE INVENTION 

GPS involves a constellation of 24 satellites placed in 
orbit about the earth by the United States Department of 25 
Defense. Each satellite continuously broadcasts a GPS sig- 
nal. This GPS signal contains an L-band carrier component 
(LI) transmitted at a frequency of 1 .575 GHz. The LI carrier 
component is modulated by a coarse acquisition (C/A) 
pseudo random (PRN) code component and a data compo- 30 
nent. 

The PRN code provides timing information for determin- 
ing when the GPS signal was broadcast. The data component 
provides information such as the satellite’s orbital position. 

The carrier component allows a receiver to easily acquire the 35 
GPS signal. 

Position determination using Conventional GPS is well 
known in the art. In Conventional GPS, a receiver makes 
ranging measurements between an antenna coupled to the 
receiver and each of at least four GPS satellites in view. The 40 
receiver makes these measurements from the timing infor- 
mation and the satellite orbital position information obtained 
from the PRN code and data components of each GPS signal 
received. By receiving four different GPS signals, the 
receiver can make fairly accurate position determinations. 45 

However, Conventional GPS only allows a user to deter- 
mine his actual location to within tens of meters. However, 
this is not suitable for applications such as attitude deter- 
mination for moving vehicles which requires extreme pre- 5Q 
cision. 

A more accurate version of GPS is Ordinary Differential 
GPS. Position determination using Ordinary Differential 
GPS is also well known in the art. It involves the same kind 
of ranging measurements as are made with Conventional 55 
GPS, except that a ground reference receiver at a precisely 
known location is utilized. Ideally, satellite ranging errors 
will affect the position determinations made by the user’s 
receiver in the same way as they will the position determi- 
nations made by the nearby ground receiver. Since the go 
location of the ground receiver is already known, the ground 
receiver can compare the position determination it has 
calculated with the actual known position. As a result, the 
ground receiver can accurately detect ranging errors. 

From these errors, the ground receiver can compute 65 
suitable corrections which are transmitted by data link to the 
user’s receiver. The user’s receiver can then apply the 


2 

corrections to its own ranging measurements so as to pro- 
vide accurate real time position determinations. 

However, even with Ordinary Differential GPS, the posi- 
tion determinations are only accurate to within several 
meters. Since, as indicated earlier, attitude determination 
must be extremely accurate, extending Ordinary Differential 
GPS to attitude determination is not feasible. 

An extremely accurate form of GPS is Carrier Phase 
Differential GPS. This form of GPS utilizes the 1.575 GHz 
carrier component of the GPS signal on which the PRN code 
and the data component are superimposed. 

Carrier Phase Differential GPS involves generating posi- 
tion determinations based on the measured phase differences 
at two different antennas for the carrier component of a GPS 
signal. However, this technique initially requires determin- 
ing how many integer wavelengths of the carrier component 
exist between the two antennas at a particular point in time. 
This is called integer ambiguity resolution. 

A number of approaches currently exist for integer ambi- 
guity resolution. However, all of them suffer from problems 
which render them unfit for applications requiring extremely 
precise attitude determinations. 

One approach is Integer Searching using redundant mea- 
surements. This involves receiving more than the standard 
four GPS signals in order to sort out the correct combination 
of integer ambiguities. The different combinations of integer 
candidates are systematically checked against a cost func- 
tion until an estimated correct set is found. However, for 
antenna separations of just a few meters, the checked 
combinations can number in the hundreds of millions. As a 
result, this approach has a propensity to arrive at wrong 
solutions. Furthermore, the configuration of the constella- 
tion of GPS satellites can only guarantee that four satellites 
will be in view at any given time. Therefore, any application 
requiring attitude determinations at any given time must not 
rely on redundant satellites for reliable resolution of the 
integer ambiguities. 

Another approach is Narrow Correlator Spacing. This 
technique involves using the PRN code of the GPS signal to 
bound the integer ambiguities. However, a significant 
amount of the time it can yield position determination errors 
of as much as several meters. This does not provide the kind 
of consistency which is required in attitude determinations 
for moving vehicles. 

Still another approach is Dual Frequency Wide-Laning. 
This approach also utilizes a second GPS signal broadcast 
by each satellite. This second GPS signal has an L-band 
carrier component (L2) transmitted at a frequency of 1.227 
GHz. The L2 carrier component and the LI carrier compo- 
nent are differenced so as to form a signal having an effective 
wavelength that is much longer than either of the two carrier 
components. From this signal, it is relatively easy to resolve 
the integer ambiguities. However, the L2 component is not 
available for civilian use. Although the denial of the second 
carrier component can be countermeasured with cross cor- 
relation technology, the performance of this type of tech- 
nology is unproven and very expensive to implement. 

One successful approach to integer ambiguity resolution 
is motion-based and has been utilized in static surveying 
applications. This approach involves taking a number of 
phase measurements while the user’s antenna and the ref- 
erence antenna are stationary. These phase measurements 
are made over a period of about 15 minutes. The phase 
measurements made during the slowly changing geometry 
of the GPS satellites will reveal the integer ambiguities. But, 
in many applications requiring attitude determination, it 
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would be impractical to require the user’s antennas to 
remain stationary for 15 minutes while the integer ambigu- 
ities are resolved. 

Several motion-based approaches for integer ambiguity 
resolution have specifically been proposed for attitude deter- 5 
mination. They involve placing antennas at various points 
around a spacecraft or airplane. The integer ambiguities can 
be resolved with rotation of the space craft or aircraft and 
taking several phase measurements. 

However, one of these methods is restricted to rotation 10 
about multiple coordinate axes. This approach is not prac- 
tical since, for example, simple rotation of an aircraft on the 
runway prior to takeoff is not sufficient to resolve the integer 
ambiguities in this method. 

Another one of these methods is limited to small angle 
rotation of the moving vehicle with respect to the GPS 
satellite geometry. However, large angle rotation provides 
significantly more information for resolving the integer 
ambiguities. Since the small angle rotation method is not 20 
able to take advantage of such information, it does not 
provide as accurate a resolution of the integer ambiguities as 
would an approach utilizing large angle rotation. 

In still another of these methods, the integer ambiguities 
can only be resolved if the airplane or spacecraft has a 25 
constant rate of rotation. This method is therefore unneces- 
sarily limiting. 

SUMMARY OF THE INVENTION 

30 

The foregoing and other problems are solved by a GPS 
attitude receiver in accordance with the present invention 
which determines the attitude of a moving vehicle. The GPS 
attitude receiver is used in conjunction with a first, a second, 
a third, and a fourth antenna mounted to the moving vehicle. 35 
Each of the antennas receives a plurality of GPS signals that 
each include a carrier component. For each of the carrier 
components of the received GPS signals there is an integer 
ambiguity associated with the first and fourth antennas, an 
integer ambiguity associated with second and fourth anten- ^ 
nas, and an integer ambiguity associated with the third and 
fourth antennas. 

The GPS attitude receiver measures phase values for the 
earner components of the GPS signals received from each of 
the antennas at a plurality of measurement epochs during an 45 
initialization period and at a measurement epoch after the 
initialization period. In response to the phase values mea- 
sured at the measurement epochs during the initialization 
period, it computes baseline vectors representing positions 
for the first, second, and third antennas with respect to the 50 
fourth antenna at an initial one of the measurement epochs 
during the initialization period. 

In response to the computed baseline vectors and the 
phase values measured at the measurement epochs during 
the initialization period, the GPS attitude receiver computes 55 
integer ambiguity resolution values representing resolution 
of the integer ambiguities. Then, in response to the com- 
puted integer ambiguity resolution values and the phase 
values measured at the measurement epoch after the initial- 
ization period, it computes values defining the attitude of the 60 
moving vehicle at the that measurement epoch. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects of the invention will 65 
become more apparent on reading the following detailed 
description and upon reference to the drawings, in which: 
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FIG. 1 shows a general view of a GPS attitude system in 
accordance with the present invention for determining the 
attitude of a moving vehicle; 

FIG. 2 provides an illustration of how integer ambiguities 
arise at an initial epoch which are then resolved during an 
initialization period required for generating the GPS attitude 
determinations for the moving vehicle; 

FIG. 3 provides an illustration of an epoch after the initial 
epoch of FIG. 2; 

FIG. 4 provides an illustration of the vector relationships 
associated with the integer ambiguities; 

FIG. 5 shows rotation of the attitude antennas of the 
attitude system of FIG. 1 about a single axis of the local 
horizontal coordinate system during the initialization period 
required for making the GPS attitude determinations; 

FIG. 6 shows the vectors representing the surveyed posi- 
tions of the antennas which are mounted on the moving 
vehicle with respect to the body coordinate system of the 
moving vehicle; 

FIG. 7 shows the rotation of the body coordinate system 
of the moving vehicle with respect to the local horizontal 
coordinate system; 

FIG. 8 provides a detailed illustration of the GPS attitude 
system of FIG. 1; 

FIG. 9 shows the method employed by the GPS attitude 
system of FIG. 8 for resolving the integer ambiguities; 

FIG. 10 shows the method employed by the GPS attitude 
system of FIG. 8 for making the GPS attitude determinations 
once the integer ambiguities are resolved. 

DETAILED DESCRIPTION OF THE 
INVENTION 

FIGS. 1-10 provide illustrations of the invention 
described herein. In these figures, like components are 
designated by like numerals. 

General Description of System and Method 

FIG. 1 shows an airplane (or spacecraft) 21. Multiple GPS 
satellites 24(1 )-(/') at known orbital positions are in view and 
broadcast GPS signals 25(1 )-(/). 

Mounted on airplane (or spacecraft) 21 is GPS mobile 
attitude system 37 which includes GPS attitude receiver 33, 
GPS top side antenna 34, and GPS attitude antennas 
35(l)-(3). Each of the components 33, 34, and 35(l)-(3) of 
the GPS attitude system 37 is mounted to the airplane (or 
spacecraft) 21 and is therefore mobile. Each of the antennas 
34 and 35(l)-(3) receives GPS signals 25(1 )-(/) and is 
coupled to attitude receiver 33. 

The GPS signals 25(l)-(/‘) are LI C/A code GPS signals. 
In other words, they contain an LI carrier component, a C/A 
PRN code, and a data component. 

The LI carrier signal is a sinusoidal wave transmitted at 
a frequency of 1.575 GHz. In the preferred embodiment, the 
LI carrier signal allows the GPS attitude receiver 33 to 
easily acquire the GPS signals 25(1 )-(/). And, as is dis- 
cussed later, it also allows the GPS attitude receiver 33 to 
compute precise position determinations for airplane (or 
spacecraft) 21 using Carrier Phase Differential GPS. 

The PRN code provides timing information enabling GPS 
attitude receiver 33 to make Conventional GPS position 
determinations for use in making the later discussed attitude 
determinations. The PRN code signal comprises a series of 
variable width pulses broadcast at a frequency of 1.023 
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MHz. Each of the GPS satellites 24(l)-(/) transmits its own 
unique PRN code which enables GPS attitude receiver 33 to 
easily identify and separate the various GPS signals it 
receives. 

The GPS attitude receiver 33 generates internally the 5 
same PRN codes at substantially the same time as do GPS 
satellites 24(l)-(/). The GPS attitude receiver 33 compares 
the PRN codes that it generates with the corresponding PRN 
codes received from the GPS satellites 24(1 )-(/). The phase 
difference needed to match the received and generated PRN 10 
codes is then computed in terms of time. 

The computed phase difference represents the time- it 
takes for the PRN codes of the GPS satellites 24(1 )-(/) to 
travel to the antennas 34 and 35(l)-(3). From the measured 
phase differences, the ranges to the broadcasting GPS sat- 15 
ellite 24(1 )-(/) can be established. With four such ranging 
measurements and the information contained in the data 
components of the GPS signals 25(1 )-(/), the positions of 
antennas 34 and 35(l)-(3) can be computed. 

The information in the data component of each of the GPS 2Q 
signals 25(1 )-(/) includes the orbital position of whichever 
of the GPS satellites 24(1 )-(/') has broadcast it. This infor- 
mation is provided as a bit stream with a frequency of 50 bits 
per second. 

The initialization of GPS attitude receiver 33 for making 25 
attitude determinations involves integer ambiguity resolu- 
tion. Integer ambiguity resolution is the process of deter- 
mining, at a particular point in time, the number of integer 
wavelengths of the carrier component of a GPS signal 
25(1 )-(/') which lie between a given pair of antennas in the 
direction of the broadcasting GPS satellites 24(l)-(/). 30 

FIG. 2 shows the situation for the first measurement epoch 
of the initialization period. It provides an illustration of how 
the later discussed integer ambiguities n 25(l)/34/35(m) arise. 

Specifically, GPS satellite 24(0 (he. the \ lh of the GPS 
satellites 24 ( 1 )— (/")) broadcasts the carrier component of GPS 
signal 25(0 (he. the \ th of the GPS signals 25(l)-(/)) in the 
direction of antennas 34 and 35(l)-(3). The integer ambi- 
guity n 25(0/34/35(m) of GPS signal 25(0 is associated with top 
side antenna 34 and the m** of the antennas 35(l)-(3). 

40 

Receiver 33 is configured to make phase measurements 
for the acquired GPS signals 25(1)— (/). Each measurement 
includes both a fractional wavelength phase component <£ /r 
and an integer wavelength phase change component <b int . 

The integer wavelength change in phase for each raw 45 
phase measurement is kept track of by receiver 33 as of the 
time the GPS signals 25(1 )-(/) was first acquired. In the 
preferred embodiment, the phase measurements are made by 
receiver 33 at a rate in the range of 1-10 Hz with each time 
step or cycle being a measurement epoch. 50 

Receiver 33 generates internally its own carrier compo- 
nent for phase comparisons with the received carrier com- 
ponent of GPS signals 25(1 )-(/). These carrier components 
are biased because at each measurement epoch the receiver 
33 has clock synchronization error AT 33 and the signal 55 
generator of GPS satellite 24(0 has a clock synchronization 
error AT 24(0 . 

As shown in FIG. 2, the unknown range r 24(0/34 between 
the transmit antenna of GPS satellite 24(0 and antenna 34, 
at the initial epoch of the initialization, includes the phase 60 
component <£ 25 ( 0/34 measured by receiver 33 and the 
unknown integer component n 25(f)/34 of GPS signal 25(0- 
And, the unknown range r 24( . )/35(m) at the initial epoch 
between the transmit antenna of GPS satellite 24(0 and the 
m ih of the antennas 35(l)-(3) includes the phase component 65 
O 25(0/35(m) measured by receiver 33 and the unknown inte- 
ger component n 25(0/35(m) of GPS signal 25(0- 


The unknown integer components n 25(0/34 and n 25(0/35(m) 
which are assigned at the initial epoch remain constant 
throughout the initialization process and the subsequent 
attitude determination process. This is illustrated in FIG. 3. 

FIG. 3 shows an epoch after the initial epoch. This second 
epoch could be during or after the initialization period. Each 
of the measurements $ 25( *)/34 ^d $ 25 ( 0 / 3 5 (m) will have 
changed since the initial epoch. This is due to the fact that 
the fractional component <£ /r and integer wavelength change 
component O int which make up the identified phase mea- 
surements have changed since the initial epoch. However, 
the assigned integer components n 25(0/34 and n 25(0/35(m) 
have not changed. 

The relationship between $ 25 ( 0/34 ^d n 25(0/34 and the 
relationship between $ 25( - )/40 and n 25(0/40 are provided as 
follows in Equations (1), and (2) respectively: 


^ > 25(iy34- r 24(iV34“ n 25(0/34+^33"^24(i) U) 

®25(iy35(my =r 24(i)/35(m)— n 25(i)/3S(m) + AT 33 ~AT 2 4( m ) ( 2 ) 

Differencing Equations (1) and (2) provides the single 
difference phase relationship given as follows in Equation 
(3): 


®25( l y34/35(m) =< ^25(|-y34~ < t , 25(i)/35(m)-^ r 24(fy34/35(m) -n 25 (i')/34/35(m) (3) 

where (A) n 25a)/34/35(m) represents the unknown integer 
ambiguity for GPS signal 25(0 associated with antennas 34 
and 35(m) and (B) Ar 24(/)/34/35(m) represents the difference in 
the unknown ranges r 24(i)/34 and r 24(i)/35(m) . 

In order to resolve the integers ambiguities n 25(l - )/34/35(m) 
properly. Equation (3) must be manipulated so as to include 
the baseline vector relationships which are associated with 
the ranges r 24(l)/34 and r 24(0/35(m) . These relationships are 
shown in FIG. 4 and are established with respect to the 
runway or local horizontal orbit coordinate system 46 which 
is defined by the along track AT, cross track CT, and altitude 
AL coordinate axes. 

The baseline vectors b 35(1) , b 35(2) , and b 35(3) respectively 
represent the unknown positions of attitude antennas 
35(l)-(3) with respect to antenna 34 in the local horizontal 
coordinate system 46 at the initial epoch of the initialization 
period. The baseline vector b 35(m) [3x1] (i.e. the m th of the 
baseline vectors b 35(1) , b 35(2) , and b 35(3) ) is provided as 
follows in Equation (4): 


&35(«) = 


btfimVAT 
blftmyCT 
[ &35(m)/AL 


(4) 


where b 35(m)/A:r , b 35(m)/CT , and b 35(m)/AL are respectively the 
unknown distances between antennas 35(m) and 34 in the 
along track AT, cross track CT, and altitude AL directions. 

The direction to GPS satellite 24(0 (he. the i th of GPS 
satellites 24(1 -(f)) in relation to antenna 34 is represented by 
the known unit direction vector s 24(0/34 [3x1] provided as 
follows in Equation (5): 


$ 24 ( 0 / 34/47 


*524(0/34 = 


*S’24(0/34/CT 


L $24(0/34/4L 


(5) 


where s 24(l)/34/Ar , s 24(/)/34/cr , s 24(i)/34/AL are respectively the 
unit distances to GPS satellite 24 (0 in the along track AT, 
cross track CT, and altitude AL directions. This vector is 
computed by receiver 33 for a GPS satellite 24(0 from the 
satellite position information contained in the data compo- 
nent of the associated GPS signal 25(0 and from the coarse 
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position fix generated by receiver 33 for antenna 34 using 
Conventional GPS positioning. 

From the preceding vector relationships shown in FIG. 4 
and developed in Equations (4) and (5), the following 
mathematical relationship is provided in Equation (6): 


Af 24(«-)/34/35(m)-^24<0 7 ^35(m) 


( 6 ) 


Combining Equation (6) with Equation (3) results in the 
following relationship in Equation (7): 


^ > 25(«)/34/35(m) - ^24(i) T ^3SCm)~ W 2SC*y34/35(m) 


(7) 
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The integer ambiguities n 25(0/34/35(m) can be computed 
during an initialization period using two different 
approaches. The first approach requires that the airplane (or 15 
spacecraft) 21 remain stationary during the initialization 
process. The second is motion-based. 

The static method is similar to that used in surveying 
applications. After several epochs of measuring <t> 2 5 (o /34 
0 2 5 (z .) /35 ( m ), receiver 33 can generate equations from Equa- 
tion (6) which are stacked in matrix form for solving the 20 
integer ambiguities n 25(0/34/35(m) . 

In order to ensure greater accuracy for the computed 
values n 25(l)/34/35(m) , receiver 33 employs more than the 
minimum number of epochs needed to compute these val- 
ues. As a result, the system benefits from the over-deter- 25 
mined set of unknowns. 

Furthermore, receiver 33 makes measurement epochs 
over a large enough time period to allow the slowly chang- 
ing GPS satellite geometry to reveal the integer ambiguities 
n 25 ( 0 / 34 / 3 s (m) * This typically requires approximately fifteen 30 
minutes. 

Additionally, where possible, phase measurements <h 25(0/ 

34 and ^ 25 ( j)/ 35 ( m ) for additional GPS signals 25(i) are made 
by receiver 33. Again, the system benefits from the over- 
determined set of unknowns. 35 

The second approach to resolving the integer ambiguities 
n 25 (o/ 34 / 35 (m) requires large angle rotation of the antennas 
35(l)-(3) about at least one of the axis of the local horizontal 
coordinate system 46, as described in Cohen and Parkinson, 
“Aircraft Applications of GPS-Based Attitude Determina- 40 
tion: Test Flights on a Piper Dakota”, ION GPS-92, Albu- 
querque, N. Mex., September 1992, which is hereby 
expressly incorporated by reference herein. FIG. 5 shows the 
vector relationships for such a rotation. 

In FIG. 5, antennas 35(l)-(3) rotate about the altitude AL 45 
axis. This may take place in space, in the air, or the ground 
at runway 23. The baseline vectors b 35(1) , b 35(2) , and b 35(3) 
are unknown at the initial epoch of the initialization process. 
The vectors Ab 35(1) , Ab 35(2) , and Ab 35(3) respectively repre- 
sent the change in positions of the antennas 35(l)-(3) at a 50 
second epoch with respect to the initial baseline vectors 
^ 35 (i )> ^ 35 ( 2 )’ b 35 ( 3 ). The vector Ab 35(m) e the m th of the 

vectors Ab 35(1) , Ab 3S(2) , and Ab 35(3) ) is provided as follows 
in Equation (8): 


A6 35(m ) = 


' ^5(m)IAT 
A&35 (m)/CT 
. Abz5(m)!AL 


( 8 ) 


where Ab 35(m)M7 , Ab 35(m)/cr , and Ab 35(m)Mi are respectively 60 
the change in position of the antenna 35(m) at the second 
epoch in the along track AT, cross track CT, and altitude AL 
directions. 

The equations generated from Equation (8) at the initial 
and the second epoch can be subtracted to establish the 65 
following relationship in Equation (9): 


8 

^^ > 25(iy34/35(w) = ^24(0 7 ^35(m) (9) 

where A<J> 25(0/34/35(m) represents the change in ® 25(0 / 34 / 35 (m) 
between the initial epoch and the second epoch. 

The equations generated from Equation (9) may be 
stacked at a number of epochs after the initial epoch to solve 
for the vectors Ab 35(m) . Thus, the vectors Ab 35(m) may be 
simultaneously computed at each of these epochs without 
resolving the integer ambiguities n 25(i)/34/35(m) . 

The antennas 35(l)-(3) are fixed to the airplane (or 
spacecraft) 21. Thus, the following constraint relationship 
may be imposed on the baseline vectors b 35Cy) and b 35(z) (i.e. 
the y th and z th of the vectors b 35(1) -b 35(3) ) as follows in 
Equation (10): 

(^35(y)+4&3s ( y J )(fo35 (z )+Afe 3 5 (z) )=fe 3 5 ( y) 7 lb3 5(z) (10) 

However Equation (10) can also be mathematically 
expressed as follows in equation (11): 

(^35(y) + ^^350»))(^35(2)' S ‘4fc 3 s (z) )=& 3 5 t> , ) 7 t 3 s (z j+Afe 35(z) :r fe 3 5 (y) + 

4^35 (y) ^35(20+^3500 ^35(z) (H) 

Thus, the Equations (10) and (11) can be combined to form 
the following relationship in Equation (12): 


^ , 35(z) 7 ^350')+^ , 350') 7 ^35(z)““^ , 350') 7 A^35(z) (12) 

Equation (12) can be stacked by receiver 33 in matrix 
form to provide equations at each epoch after the initial 
epoch for solving the unknown vectors b 35(w) . This includes 
the situations where y^z and where y=z. 

For greater accuracy more than the minimum number of 
epochs needed to calculate the baseline vectors b 35(/n) should 
be employed by receiver 33. As a result, receiver 33 can 
generate additional equations from Equation (12) for simul- 
taneously solving the over-determined set of unknown base- 
line values. 

Once these baseline values are computed, receiver 33 then 
computes the integer ambiguities. This is done in the fol- 
lowing manner. 

The issue of the flexibility of airplane (or spacecraft) 21 
must be first taken into account for computing the integer 
ambiguities n 25(0/34/35(m) . Specifically, the difference in rela- 
tive position of antennas 35(l)-(3) with respect to antenna 
34 due to wing flexure can amount to a significant fraction 
of a wavelength. It is important to relate the flexure of 
airplane (or spacecraft) 21 to (A) the known position vectors 
^ 0 / 35 ( 1 )’ ^o/ 35 ( 2 ) > and ko/ 35 ( 3 > for the antennas 35(l)-(3) in the 
absence of flexure, (B) the integer ambiguities n 25(i)/34/35(m) , 
and (C) the phase values <F 25(0/34 and <& 25(0/35(m) measured 
for each of the GPS signals at each epoch. 

FIG. 6 shows the vector relationships associated with 
antennas 35(l)-(3) in a body coordinate system 47 defined 
by the airplane (or spacecraft) 21 and the coordinate axes X, 
Y, and Z. The known position vectors ko,^^, k^^, and 
^ 0 / 35 ( 3 ) * n body coordinate system 47 are precisely 
surveyed in the absence of flexure. The known vector 
ko/ 35 (m) [3x1] (i.e. the m th of the known vectors ko^^, 
ko/ 35 ( 2 )> and ifo /35(3 )) is provided as follows in Equation (13): 


%35(m) = 


£o/35(m)/X 
*0/3 5(m)IY 
*0/35 (m)/Z 


(13) 


where ko /35(m)/x , ^ K5imVY , and ko /35(m)/z are respectively the 
known distances between antennas 35(m) and 34 in the X, 
Y, and Z directions. 

The relationship of each known position vector k 0/35(w) to 
flexure of airplane (or spacecraft) 21 is given as follows in 
Equation (14): 



9 
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^5(m)-^3S(m)+^ex/35(m)S/ (14) 

where (A) k 35(m) is the estimate for the position of antenna 
35(m) during flexure of airplane (or spacecraft) 21, (B) 
V*» 5 (m) is the m"‘ of the flexure vectors k flexn5m , k^ 5(2) , 
and kji ex/3 5(3) and is associated with antenna 3 5(m), and (C) 

8f is the state of flexure of airplane (or spacecraft) 21 during 
the initialization process. 

The flexure vectors k^ /35(1) , k^ x/33(2) , and kj, ex/35m 10 
respectively associated with antennas 35(l)-(3) are given as 
follows in Equations (15)— (17): 


kjlexJ 35(1) = 


0 

0 

-1 


kjlexJl 5(2) = 


0 

0 

-1 


(15) 

15 

(16) 

20 


kjlejl 35(3) = 


0 

0 

0 


(17) 


FIG. 7 shows the vector relationships associated with 
antennas 35(1)— (3) as the body coordinate system 47 rotates 
about the local horizontal coordinate system 46. The 
unknown vectors x 35(1) , x 35(2) , and x 35(3) respectively rep- 
resent the unknown positions of attitude antennas 35(1), 
35(2), and 35(3) from antenna 34 with respect to the local 
horizontal coordinate system 46. The unknown vector x 35(m) 
[3x1] (i.e. the m th of the unknown vectors x 35(1) , x 35(2) , and 
X 35 ( 3 >) is provided as follows in Equation (18): 


25 


30 


*35 (m)fAT 


*3S(m) = 


*35(m)/CT 


L *35(m)/AL 


(18) 


35 


where x 35(m)/A7> x 35(m)/C7> and x 35(m)/Ai are respectively the 
unknown distances between antennas 35(m) and 34 in the 40 
along track AT, cross track CT, and altitude AL directions. 

From the preceding vector relationships shown in FIG. 7 
and developed in Equations (5) and (18), the following 
relationship is provided in Equation (19): 


^ r 24( ; y34/35(m)^ 7 24(f>*35(m) 


(19) 


Combining Equations (19) and (3) results in the following 
relationship in Equation (20): 


^25(0/34/35(m)^24(0 7 *3S(»0 n 25(i)/34/35(m) 


50 

( 20 ) 


Furthermore, the attitude matrix A represents the estimate 
of the attitude of the airplane (or spacecraft) 21 at a 
particular epoch. The matrix A is established from the 
rotation of the body coordinate system 47 of airplane (or 55 
spacecraft) 21 with respect to the local horizontal coordinate 
system 46 at the epoch. The transpose of this matrix is 
provided as follows in Equation (21): 


Ax/at 


A T = 


Ax/a' 

Axial 


An at 
Ancr 
Ay/al 


Ayjat 

Ancr 

A-z/al 


( 21 ) 60 


where each element of the matrix represents the rotation of 
a coordinate of the body coordinate system 47 with respect 65 
to a coordinate of the local horizontal coordinate system 46 
at the measurement epoch. 


10 

The following relationship may then be established in 
Equation (22): 

*35(m)“^ r ^35(m) (22) 

As a result, the vector A r k 35(m) represents the position of 
antenna 35(m) relative to antenna 34 in the local horizontal 
coordinate system 46. Therefore, combining equation (19) 
with Equation (20) results in the following relationship 
provided by Equation (23): 


< t ) 25(iy34/35(m)-^24(0 r ^ 7 ^35(m) -n 25(0/34/35(m) (23 ) 

Furthermore, the attitude matrix A may be expressed as 
follows in Equation (22): 

A=&AA 0 (22) 

where (A) 5A [3x3] is an attitude correction matrix of small 
angle rotations, and (B) A 0 is the estimate for A from the 
previous iteration. The correction matrix 5A is expressed as 
follows in Equation (23): 


5A=/+0* (23) 

where (A) I [3x3] is an identity matrix, and (B) 0* [3x3] is 
the skew symmetric matrix at each measurement epoch 
associated with the unknown vector 80 of small angle 
rotations at the same epoch. 

The unknown vector 80 [3x1] can be expressed as 
follows in Equation (24): 



60z 


(24) 


where 8©*, 80 y , and 80 2 respectively represent the 
unknown small angle rotations about the X, Y, and Z 
coordinates of the body coordinate system 47. 

The skew symmetric matrix 0* associated with the vector 
8© can be expressed as follows in Equation (25): 


0* = 


ro 

56 z 
— 58 k 


-5 Q z 
0 

50x 


66y 

-50x 

0 


(25) 


Combining Equations (14), (22), (23), with Equation (21) 
results in Equation (26): 


kari5(m) T &A Q $24(r > +kji e ^ 5{m ) T AQ$24W§f+kjlex/3S(m) T & ^ 0 ^ 24 ( 0 ^/“ 

n 25(i)/34/35(m)-^25(i)/34/35 4 o?24(0 (26) 


However, the product of 0* and 8f is assumed small and 
therefore the term k //ejc/35(m) r 0 x A o s 24(O 8f can be neglected 
in Equation (26). Furthermore, the dot product of the matrix 
0* and the vector ko /35(m) can be expressed as the dot 
product of the matrix K 0/35(m) * and the vector 80 where the 
matrix K 0/35(m) * is represented as follows in Equation (27): 


ro 


-*0/35(m)Mr *0/35 (m)/AT 


(27) 


*0/35(m) “ 


*0/35(m)/CT 0 


“*0/35(m)/CT 


L ~*0/35(m)M *0/35(m)/A 0 J 

Thus, Equation (26) may be rewritten as follows in Equation 
(28): 


^24(0 7 ^O^CV35(m) J: ^|*/T^5(m) 7 4(>?24(0^/~ n 25(0/34/35Cm) := ^ ) 25(iy34/ 
35(m)“*o/35( m ) 2 A (> ?24(0 (28) 

Note that Equation (28) does not require or assume that 
airplane (or spacecraft) 21 is rotating at a constant rate of 
rotation during the initialization period. Thus, this equation 
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may be applied to the situation when airplane (or spacecraft) 
21 is rotating at a constant rate of rotation and to the situation 
when it is not rotating at a constant rate of rotation. 

Equation (28) can then be stacked over the N measure- 
ment epochs of the initialization period so as to from the 
following set of simultaneous Equations defined by the 
matrix system (29): 


1 

■ i 


(29) 

0 0 1 f) 1 /' 


fie<» 


‘ 8 Vd) ' 

o : i : i i 


68(4) 



: m) \ i /?> i i 



= 

Svw 

: 0 1 : l : 


80WJ 



,0 0 HW 1 /*> 1 /. 


n 

i ■ 




where (A) each small angle rotation vector 80 (<7) is the small 
angle rotation vector at the q th of the N epochs, (B) the 


12 

of the antennas 35(l)-(3), each of the matrices H c<?) will be 
18x3. Combining together perhaps 60 phase measurements 
(i.e. N=60) O 25(0/34 and O 25(0/35C m) over a thirty second 
span, the total size of the system is 1080x199. The resulting 
5 storage is approximately 1 MB of RAM. 

However, the sparsity of the coefficient matrix for the left 
side of the set of Equations (29) can be taken advantage of 
in solving for the vector of unknowns 80 (<?) , 5f and n with 
10 a minimum of storage space. This can be done with block 
Cholesky decomposition as described in Golub and van 
Loan, “Matrix Computations”, 2nd Ed., Johns Hopkins 
University Press, Baltimore, Md., 1989, pp. 141-145, which 
is hereby expressly incorporated by reference herein. 

15 An especially efficient method of implementing the block 
Cholesky decomposition is provided next. This first requires 
re-constructing the set of Equations defined by the matrix 
system (29) as follows in the symmetric matrix system (33): 


Ftfwwi) o 


o I tf» T 


(33) 


L^- 

■ m 


mm) 

o 


_ o m T m 
m ) ■ ■ ■ 7/W) 


mv 


m T 

1 




V ij 1 

60(4) 



50W 


y^) 

■ n ■ 


. Vn . 


vector n comprises the integer ambiguities n 25( . )/34/35(m) , and 
(C) the I matrices are identity matrices. Furthermore, the 
vector 8\j/ (9) in Equation (29) is the vector of differential 
phase residuals at the q th of the N measurement epochs and 40 
is given by Equation (30): 

= [®2?(t)/34/35(m) “ ^35(m)^(f ) ‘ s 24(r)] ( 30 ) 

where each differential phase measurement $ 25 (o/ 34 / 35 (m) C<?) 
is the differential phase at the q lh of the N measurement 45 
epochs. The matrix H (9) (i.e. the q th of the matrices H (1) 
through H (A °) in Equation (29) is the sensitivity matrix of 
changes in measured differential phase with 

respect to the small angle rotations 80 (?) at the q"* of the N 
measurement epochs and is given by Equation (31): M) 

< 31 > 

where the matrix A 0 (<?) is the attitude matrix at the q th of the 
N measurement epochs of the initialization period. The 55 
vector f^ in Equation (29) is the flexure sensitivity vector 
of changes in differential phase d> 25(0/34 , 35(m) (<7) with respect 
to the flexure state value Sf at the q th of the N measurement 
epochs and is given by Equation (32): 

60 

/ 9) = [^x/35 (m /0^24(0] < 32 > 

The vector of unknowns S0 (£?) , Sf, and n can be iteratively 
solved using the matrix system (29). However, this matrix 
system may be overwhelming for the computing power of 65 
GPS attitude receiver 33. For example, considering a sce- 
nario with 6 of the GPS signals 25(l)-(/) received by each 


where (A) the q th of the matrices H (1) through H (A ° is given 
as follows in Equation (34): 


m)= 


[ 


fq)rm ) 


] 


(B) the matrix I is given as follows in Equation (35): 


(34) 


1 = 


£ fjftTfa) 

N 

£ /4> 
4=1 


N 

£ fq)? 

4=1 

NI 


(35) 


(C) the q ih of the vectors \j / (1) through \j/ 1} is given as 
follows in Equation (36): 




(36) 


(D) the vector is given as follows in Equation (37): 


Yn = 


4=1 


£ 

4=1 


and (E) the vector n is given as follows in Equation (38): 

"=[!■] (38) 
For convenience, the matrix system (33) can be described 
in short hand notation in the following matrix system (39): 


Hx=\\f s 


(39) 
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where (A) the coefficient matrix H on the lefthand side is the 
large symmetric sensitivity matrix, (B) the column vector \|f 
on the righthand side is the large differential phase residual 
vector, and (C) the solution vector x is the large vector of 
unknowns. 

Receiver 33 performs block Cholesky factorization of the 
symmetric matrix H of the matrix system (39) so as to define 
the following symmetric matrix system (40): 

LL r x=\\f (40) 


where LL r equals the matrix H. The symmetric matrix L is 
given as follows in Equation (41): 


L = 


lP> 0 


. 0 I 0 


lW 


(41) 


15 


20 


25 



. . o m 

GW • ■ ■ GW 



30 


where (A) L (<?) is the q th of the individual matrices L (1) 
through L (a ° and is given as follows in Equation (42): 


lw s M ifiw'iw 

35 

(B) the o[ h of the matrices G (<7) through G (A ° is given as 
follows in Equation (43): 


Q(g)T =L (q)-l^q)T 


(43) 

40 


and (C) the matrix L rt is given as follows in Equation (44): 



. N 

/- Z GWGWT 
q=l 


(44) 


45 


Then receiver 33 utilizes Cholesky block back substitu- 
tion for computing the vector x in matrix system (40). This 
is done by receiver 33 by first solving for the vector y given 
in the following relationship (45): 


y=L"V 


(45) 


where (A) y=L r x. The vector y is expressed as follows in 
Equation (46): 


y(l) 


55 

(46) 


yW 


60 



where (A) y (<?) is given as follows in Equation (47): 65 


14 

and (B) y n is given as follows in Equation (48): 

, / _ N _ \ (48) 

J 

Once the vector y in the matrix system (45) has been 
computed by receiver 33, then the vector x of unknowns can 
be solved for with the following relationship in Equation 
(49): 


x=L~ T y (49) 

where (A) the vector 80 (9) included in the vector x is given 

as follows in Equation (50): 

WW=&>- T (yW^rn) (50) 

(B) the vector n included in the vector x is given as follows 
in Equation (51): 

n=L~ T y n (51) 

Then, the vectors 80 (9) of small angle rotations and the 

integer ambiguities n 25(0/34/35(m) can be obtained through 
Equation (38). Thus, through block Cholesky decomposi- 
tion, the vectors S0 (<?) of small angle rotations, the integer 
ambiguities n 25(0/34/35(m) , and the flexure state value 8f can 
all be solved for by receiver 33. 

However, in order to iteratively compute the vector n, an 
initial estimate for the attitude matrix A 0 (<?) at each of the N 
measurement epochs of the initialization period must be 
computed so that initial estimates for the matrices H (9) and 
the vectors and \{/ 9) can be determined. This is done in 
the manner described in Cohen, Cobb, and Parkinson, “Two 
studies of High Performance Attitude Determination Using 
GPS: Generalizing Wahba’s Problem for High Output Rates 
and Evaluation of Static Accuracy Using Theodolite”, ION 
GPS-92, September 1992, and in Markley, “Attitude Deter- 
mination using Vector Observations and the Singular Value 
Decomposition”, Journal of the Astronautical Sciences, Vol. 
36, No. 3, July-September 1988, pp. 245-258, both of 
which are hereby expressly incorporated by reference 
herein. 

Receiver 33 does this by first computing the positions of 
the antennas 35(l)-(3) with respect to antenna 34 in the local 
horizontal coordinate system 46 at each of the N measure- 
ment epochs of the initialization period. Thus, the position 
for the m** of the antennas 35(l)-(3) at the q lh measurement 
epoch can be given by the relationship in Equation (52): 

AW* (52) 

where (A) Ab 35(m) (<?) is the q th position change vector for m** 
of the antennas 35(l)-(3) computed from the stacked set of 
equations generated from Equation (9) at each epoch of the 
N epochs as described earlier, and (B) b 0/35(m) is computed 
from the stacked set of equations generated from Equation 
(12) as described earlier 

Then receiver 33 can compute an intermediate matrix at 
each of the N measurement epochs of the initialization 
period so that the intermediate matrix D (9) for the q th of the 
measurement epochs is described as follows in Equation 
(53): 


D^=KB iq)T (53) 

where (A) the matrix K is given as follows in Equation (54) 

^=[^ 35 ( 1 )^ 35 ( 2 )^ 35 ( 3 )] ( 54 ) 

and (B) the matrix B (9) is given as follows in Equation (55): 


/*)=L(*)-y9) 


(47) 
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(55) 

Receiver 33 then performs the singular value decompo- 
sition (SVD) of each of the matrices D (<?) to generate the 
orthonormal matrices U ig) p and V (<?) 0 in accordance with 
the following relationship in Equation (56): 

£>(?) = ifgtigVg (56) 

where (A) U (9) £) and V ig) D are the orthonormal matrices at 
the q th of the N epochs, and (B) 10 

2$ 


/= Z . Z W35( m )/24(0( Ar 24(i)/34/35(m) “ 

m=l (=1 

(^ 0 / 35 (m) + fyfe*/35(m)(S/o + 85/))fiAAoS24(i)) 2 


where (A) the estimate A 0 is the attitude matrix A computed 
at the previous measurement epoch, (B) the estimate Sf 0 is 
the flexure state value 8f computed for the previous epoch, 
and (C) S8f is the flexure state correction value computed at 
each epoch. The relationship of 5f, 5f 0 , and 8Sf is given as 
follows in Equation (60): 


is a diagonal matrix of positive singular values at the q th of 15 
the N epochs. 

The signs of the third column (associated with the small- 
est singular value) of both orthonormal matrices U (?) G and 
V ((?) G are adjusted by receiver 33 to yield the proper matrices 
U (g \ and W ig \. The optimal attitude matrix A (9) 0 at the q lh 
measurement epoch is then computed by receiver 33 from 20 
the following relationship in Equation (57): 

(57) 

Once the very first initial guess for the attitude matrix 25 
A (<?) 0 is computed, the vector n is iteratively computed by 
taking the vectors of small angle rotation 80 ((?) computed 
from the current iteration and substituting them into the 
Equation 25 to form the skew matrices ©* (9) . The computed 
matrices 0* (9) are then substituted into the Equation (23) to 30 
form the attitude correction vectors 8A (9) . Then, the current 
estimate for the attitude matrices A 0 (<7) is computed and 
squared off in accordance with Equation (22). These esti- 
mates are then substituted back into the matrix system (33) 
for computation of the next estimate of the vector n. This is 35 
repeated until the computed values which represent the 
unknown vectors 80 (9) , the integer ambiguities n 25(0/34/ 
35(m), and the flexure state value 8f converge. Thus, the 
iteratively computed values for the integer ambiguities 
n 25(z)/34/35(m) represent resolution of the integer ambiguities. 40 
As was indicated earlier, these integer ambiguities n 25(0/34/ 

3 5(m) can be computed for when airplane (or spacecraft) 21 
is at a constant rate of rotation and for when it is not at a 
constant rate of rotation because of the way in which 
Equation (28) is derived. 45 

Once the integer ambiguities n 25(0/34/35(m) have been 
resolved, the initialization process is over and attitude solu- 
tions for airplane (or spacecraft) 21 (i.e., the attitude matrix 
A at each epoch after the initialization period) can then be 
computed. 50 

A complete attitude solution can be generated by receiver 
33 at each epoch after the initialization period by utilizing 
the differential ranges Ar 24(0/34/35(m) which can be computed 
from Equation (3). This is done by minimizing the following 
quadratic cost function provided in Equation (58): 55 

1 7 T - . (58) 

J= Z . Z W35( m )/24(0 (Ar 2 4(i)/34/35(m) “ *35(mA^4(i)) 
ro= 1 i=l 

where w 35(m)/24( i)represent the optional measurement 
weighting associated with antenna 3 5(m) and GPS satellite 60 
24(i)- Equation (53) is derived similarly to Equation (21). 

Starting with an assumed estimate A 0 for the attitude 
matrix A for the current attitude matrix A and an assumed 
estimate Sf 0 for the current flexure state value, a better 
estimate for the current attitude matrix A may be obtained by 65 
linearizing Equation (58) about the estimates A 0 and 8f 0 as 
follows in Equation (59): 


S/=5/ 0 +55/ (60) 

After combining Equations (22)-(24) and (60) with Equa- 
tion (59) (in a similar way to which Equations (22)-(24) 
were combined with Equation (21) to produce Equation 
(28)), the attitude cost function can be expressed as follows 
in Equation (61): 

1 i (61) 

/= Z Z W35( m )/24(o(8r 2 4(i)/34/35(m) “ 
m = 1 i=l 


^5(^^24(088/- 5 2 4( i yto 7 ’/^5 (OT) 50) 2 


(B) the skew matrix K 35(m) * is given as follows in Equation 
(63): 


ro 


-*35(m)Mr 


«W B 


*35(m)/CT 
L _ ^35(m)M 


0 

&35(m)M 


and (C) the vector 8r 24(0/34/35(m) 
in Equation (64): 


&35(m)Mr 

-k$5(m)fCr 


(63) 


0 J 

at the current epoch is given 


8 ^24 (i)/34/35 (m)~^*24 (i)/34/35(m) ^3 5 (m)^ 0^24(») (64) 

The linearized cost function can now be written as a matrix 
system (65): 


i 4 ii / r i 

wUHl Hrtn H"§4" 8r IL 

where (A) W 1/2 is a diagonal matrix formed by the square 
root of each weighting factor w 35(m)/24(0 and in the preferred 
embodiment is an identity matrix, (B) the matrix H is the 
sensitivity matrix at the current epoch and is given by 
Equation (66): 


H= 


524(2) ^O 7 ” Sy (i) Aq T fytex/35(m) 


By minimizing the matrix system (65), the correction 
vector 8© and the flexure state correction value S8f may be 
iteratively computed by receiver 33 after a new set of phase 
measurements O 25(0/34 and d> 25(0/35(m) have been made for 
a given epoch. As a result, the matrix 0* may be computed 
from Equation (25) and the matrix 8A may then be computed 
from Equation (23). Using the computed matrix 8A and the 
computed matrix A 0 , receiver 33 computes a more accurate 
estimate for matrix A from Equation (22). Furthermore, Sf is 
computed using Equation (60). The estimates A and 8f from 
the previous iteration are used as the current solutions A 0 
and 8f 0 for the next iteration. 
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The new estimates for A and 8f are then computed and 
compared with the estimates from the previous iteration. 
This process is continued until the estimates for A and 8f at 
the measurement epoch converges to within a desired level. 

Thus, a flexure state value 5f is computed at each epoch 5 
after the initialization period. This is different from the 
initialization process where the flexure state value 5f is 
assumed constant over the initialization period. 

The matrix A and the computed flexure state value 5f just 
computed are then used as the initial guesses A 0 and 8f 0 for 
the next epoch. Then, the iterative process described above 
is repeated to obtain the attitude matrix A and flexure state 
value 8f for this epoch. 

Another significant advantage to this approach is that the 
integer ambiguities n 2 5 (i)/ 34 / 35 (m) of an additional GPS signal 
25(0 broadcast by GPS satellite 24(0 which was not in view 15 
during the initialization period can be resolved once it does 
become in view after the initialization period. Receiver 33 
accomplishes this by measuring 0 25(/)/34 and <b 25(0/35(m) for 
the new GPS signal 25(0 at a particular epoch after the 
initialization period. At this epoch the matrix A has already 20 
been determined by receiver 33 from the other GPS signals 
25(1 )-(/) which have had their respective integer ambigu- 
ities n 25(0/34/35(m) resolved during the initialization period. 

The calculated attitude matrix A and the phase measure- 
ments O 25(0/34 and O 25(0/35(m) are plugged into Equation 2 s 
(20) so as to solve for the rounded off new integer ambiguity 
n 25 (i)/ 34 / 35 (m>* Th en > this newly computed integer ambiguity 
n 25 c 0 / 34 / 3 5 (m) together with the phase measurements for the 
newly acquired GPS signal 25(0 may be used in computing 
the matrix A in the way just described. Thus, this technique 
results in a seamless integer hand-off so that a new initial- 30 
ization period is unnecessary. 

Also, due to shadowing or masking of one of the antennas 
34 or 35(1)— (3), receiver 33 will occasionally produce an 
error in its carrier phase tracking which results in an integer 
(or cycle) slip. An integer slip is a miscount in the integer or 35 
multiple integer change of the integer wavelength portion 
<D in , between successive phase measurements O 25(0/34 or 
<& 25 (z)/ 35 (m) ® or the carrier component of a GPS signal 25(0 
at one of the antennas 34, 35(1), 35(2), or 35(3). In other 
words, the change in the integer portion is not the actual 40 
integer change in range between these successive phase 
measurements. As a result of an integer slip, the integer 
portion of all of the succeeding phase measurements 
3*25(0/34 or 3*25(0/3 5 (m) f° T the carrier component of the GPS 
signal 25(0 received by the same antenna 34, 35(1), 35(2), 45 
or 35(3) will then be in error by the integer amount of the 
integer slip. 

Since the phase measurements 0 25(z734 or <D 25(0/35(m) for 
the GPS signal 25(0 at this antenna 34, 35(1), 35(2), or 35(3) 
will be in error, the computed attitude matrix A at each epoch 50 
will then also be in error. In order to correct such a problem, 
receiver 33 will calculate a new integer ambiguity resolution 
value n 25(l)/34/35(m) which is associated with the particular 
phase measurements 0 25(i)/34 or O 25(0B5(m) offset by the 
integer slip. This new integer ambiguity resolution value 55 
n 25 (i)/ 3 4 / 3 5 (m> is computed in the same way as was described 
above for a GPS signal 25(0 that just came into view. Thus, 
receiver 33 calculates the new integer ambiguity resolution 
value n 25(l - )/34/35(m) while continuing to use the phase mea- 
surements O 25(0/34 or ff* 25 ( z)/ 35 (m) offset by the integer slip so 60 
that the new integer ambiguity resolution value n 25(0/34/35(m) 
takes into account the integer slip amount. 

Detailed Description of Attitude System and 
Method 

65 

FIGS. 8-13 provide a detailed illustration of the GPS 
attitude system 37. The functions of the components of the 
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attitude system 37, in relation to the previously described 
equations, are better understood with reference to these 
figures. 

FIG. 8 shows the attitude mobile system 37. It includes 
the GPS attitude receiver 33 and antennas 34 and 35(l)-(3). 

Antennas 34 and 35(l)-(3) receive GPS signals 25(1)— (/). 
As was indicated earlier, the positions of antennas 35(1)— (3) 
with respect to antenna 34 are respectively given by the 
vectors x 35(1) , x 35(2) , and x 35(3) in the local horizontal 
coordinate system 46 and given by the vectors k 35(1) , k 35(2) , 
and k 35(3) in the body coordinate system 47. 

GPS receiver 33 receives the GPS signals 25(1 )-(/) from 
each of the antennas 34 and 35(l)-(3). It includes a signal 
receiving block 110, a signal processing block 111, a refer- 
ence oscillator 115, a synthesizer 116, and a computer 117. 

The signal receiving block 110 comprises four signal 
receiving stages 113(l)-(4). The signal receiving stage 
113(4) is coupled to antenna 34 for receiving the GPS 
signals 25(1 )-(/'). The signal receiving stages 113(1)— (3) are 
respectively coupled to antennas 35(l)-(3) for also receiving 
the GPS signals 25(1 )-(/). Each of the signal receiving 
stages 113(l)-(4) extracts the GPS signals 25(l)-(/) it 
receives and down converts them to an intermediate fre- 
quency for signal processing by the signal processing block 

HI 

The signal processing block 111 includes 4 multi-channel 
signal processing stages 114(l)-(4). The signal processing 
stages 114(l)-(4) are respectively coupled to the signal 
receiving stages 113(l)-(4). The signal processing stages 
114(l)-(4) are also coupled to computer 117 for receiving 
signal processing control signals from the computer 117. 
Each of the signal processing stages 114(1)— (4) separates 
(i.e. demodulates) each of the down converted GPS signals 
25(1)— (/') it receives into its carrier, PRN code, and data 
components. 

Furthermore, with the signal processing control signals 
provided by the computer 117, each of the signal processing 
stages 114(1)— (4) phase locks the carrier and PRN code 
components of each of the GPS signals 25(1 )-(/) it receives 
with the carrier and PRN code signals it generates. As a 
result, the signal processing stages 114(l)-(4) provide the 
computer 117 with information for making the earlier 
described carrier phase measurements O 25(0/34 and 3> 25(0/ 
35(m> for the GPS signals 25(1)— (;) received with the anten- 
nas 34 and 35(l)-(3) and also provides signal strength 
information for each of these GPS signals. 

The computer 117 is coupled to the signal processing 
stages 114(l)-(4). It includes a central processing unit 
(CPU) 118 and a computer memory 119. 

The CPU 118 receives from the signal processing block 
111 the information for making the earlier described carrier 
phase measurements <£ 25 ( 0/34 and <£ 2 5 (o/ 35 (m) described ear- 
lier for the GPS signals 25(1 )-(/) received with the antennas 
34 and 35(1)— (3). Furthermore, the CPU also receives from 
the signal processing block 111 the demodulated data com- 
ponents of the GPS signal 25(1 )-(/). In response to this, the 
CPU outputs attitude determinations for the airplane (or 
spacecraft) 21. 

The computer memory 119 stores the signal processing 
routine 120, the carrier phase measuring routine 121, the 
directional vector computation routine 122, the static atti- 
tude initialization routine 123, the motion based attitude 
initialization routine 124, the attitude generating routine 

125, and the integrity and integer ambiguity hand-off routine 

126. The computer memory also stores data generated from 
these routines 120-126 in the data storage area 127. The 
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CPU 118 is coupled to the computer memory 119 for vectors s 24(0/34 computed by routine 122. Routine 130 

receiving the routines 120-126 and the data in the data records these values in the data storage area 127 so that 

storage area 127. Equation (12) is stacked in matrix form. As a result, the 

The signal processing routine 120 generates the signal baseline vectors b 35(m) can be simultaneously solved and 
processing control signals for controlling the carrier and 5 stored in the data storage area 127. 

PRN code phase locking operations of the signal processing At step 138, the initial attitude guess subroutine 131 first 

block 111. These control signals are output by the CPU 118 generates the N matrices B (<?) and the matrix K according to 

and received by the signal processing block 111. Equations (54) and (55) in the manner described earlier. 

The carrier phase measuring routine 121 makes the phase Specifically, subroutine 131 generates the N vectors b 35(m) (1) 

measurements d> 25( . )/34 and d> 25(0/35(m) based on the infor- 10 through b (A ° 35(m) which form the matrix B iq) according to 
mation received from the signal processing block 111. Thus, Equation (52) from (A) the baseline vectors b 35(m) computed 

the routine 121 and the signal processing block 111 make up ^ subroutine 130, and (B) the N baseline change vectors 
the carrier phase measuring component of the receiver 33. Ab ^ throu S h ^ 35(m) w computed by subroutine 130. 

As was indicated earlier, each of these carrier phase mea- Then, as shown in step 140, the subroutine 131 computes 

surements includes both a fractional wavelength phase com- 15 the N matrices D (9) according to Equation (53). In other 
ponent 0 /r and an integer wavelength phase change com- words, the N matrices D (q) are computed from (A) the N 

ponent These phase measurements are used by receiver matrices B (<7) , and (B) the matrix K. 

33 for making the Carrier Phase Differential GPS attitude In step 142, the subroutine 131 computes the N orthonor- 

determinations. As was discussed earlier, these carrier phase mal matrices U (<?) D and the N orthonormal matrices V ( ^ )r D 

measurements are made by the routine 121 at the rate of 20 according to Equation (56) from the N matrices D (<?) com- 

approximately 1-10 Hz. puted by subroutine 131 in step 140. 

The unit directional vector computation routine 122 com- Then, at step 144, the subroutine flips the signs of the third 

putes the vectors s 24(0/34 in the manner described earlier. column (associated with the smallest singular value) of the 

Thus, these vectors are computed from (A) the satellite N orthonormal matrices U C<?) D and the N orthonormal matri- 

orbital positions received in the data components of the GPS 25 ces to yield the N proper matrices U (<?) + and N proper 
signals 25(1 )-(/), and (B) the location of reference antenna matrices V (<?) + . 

34 in the coordinate system used to define the satellite orbital And finally, at step 146, the subroutine 131 computes the 

positions computed by routine 122 from Conventional GPS. n initial guess attitude matrices A (<?) 0 according to Equation 

The static attitude initialization routine 123 when selected (57) from (A) the N matrices U (<?) + , and (B) the N proper 

computes the initialization values n 25( - )/34/35(m) from Equa- 30 matrices V (q \. 

tion (7) in the manner described earlier. Thus, routine 123 is At step 148, the integer ambiguity resolution subroutine 

responsive to (A) the measurements O 25(0/34 and O 2 5 C0 / 3 5(m) 132 first computes the N sensitivity matrices H (9) , the N 

made by routine 122 over several epochs, and (B) the vectors f q \ and the N vectors 5\{/ (<?) in accordance with the 

directional vectors s 24(0/34 computed by routine 122. Since matrix system (29). Thus, the N vectors are computed 
routine 123 records these values so that Equation (43) is with Equation (30) from (A) the N phase measurements 

stacked in matrix form, the values n 25(0/34/35(m) can be ^ 2 5 (o/ 34 / 35 (m) (9) for eac b GPS signal 25(i) received by 

simultaneously solved. These initialization values n 25(0/34/ antenna 35(m), (B) the vectors k 0/35(m) stored in data storage 

3 5<m) are then stored in the data storage area 127 for use by area 127, (C) the N attitude matrices AS q \ computed by 

the attitude determination routine 125. As a built in integrity subroutine 131 for the first iteration and, as described 

check, these values are checked to see that they converge to shortly, computed by subroutine 132 for the succeeding 

integer values. iterations, (D) the unit directional vectors s 24(0/34 computed 

The motion based attitude initialization routine 124 when by routine 122. The N matrices H c<?) are computed in 

selected computes the initialization values n 25(0/34/35(m) in accordance with Equation (31) from (A) the matrices 

the manner described earlier. As shown in FIG. 8, the 45 Ko /35(m * stored in the data storage area 127, (B) the N 
initialization routine 124 includes a baseline vector compu- attitude matrices A (q) 0 computed by subroutine 131 for the 

tation subroutine 130, an initial attitude guess subroutine first iteration and, as described shortly, computed by sub- 

131, and an integer ambiguity resolution subroutine 132. routine 132 for the succeeding iterations, and (C) the unit 

The operation and relationship of these routines to one directional vectors s 24(0/34 computed by routine 122. The N 

another is shown in FIG. 9. 50 vectors f (9) are computed in accordance with Equation (32) 

At step 134, subroutine 130 first computes the vectors from (A) the vectors k^ cjc/35(m) recorded in storage area 127, 

Ab 3S(m) based on Equation (9) in the manner described (B) the N attitude matrices A (?) 0 computed by subroutine 

earlier and records them in the storage area 127. Specifically, 131 for the first iteration and, as described shortly, computed 

routine 130 records in data storage area 127 the measure- by subroutine 132 for the succeeding iterations, and (C) the 

ments 0 25(r - )/34 and 0> 25(0/35(m) made by routine 121 at an 55 unit directional vectors s 24(f)/34 computed by routine 122. 
initial epoch. Then, at a number of succeeding epochs Then, at step 150, subroutine 132 constructs the matrix 

routine 130 computes the vectors Ab 35(m) from (A) the system (33). Thus, subroutine 132 constructs (A) the N 

measurements O 25(0/34 and d> 25(0/35(m) recorded at the initial matrices H (f/) in accordance with Equation (34), (B) the 

epoch and made at these succeeding epochs, and (B) the unit matrix I in accordance with Equation (35), (C) the N vectors 

directional vectors s 24(0/34 computed by routine 122. 60 \jjr (9) in accordance with Equation (36), and (D) the vector 

Then, at step 136, baseline vector computation routine m accordance with Equation (37). 

130 computes the baseline vectors b 35(m) from Equation In step 152, subroutine 132 performs block Cholesky 

(12). Specifically, these values are generated from (A) factorization of matrix H in matrix system (39) to construct 

measurements d> 250)/34 and d> 25(0/35(m) made by routine 121 the matrix system (40). Thus, it computes (A) the N matrices 

at the epochs employed, (B) the vectors Ab 35(m) computed 65 h iq) in accordance with Equation (42), (B) the N matrices 
by subroutine 130 from Equation (9) at each of the epochs G (q)T in accordance with Equation (43), and (C) the matrix 

employed after the initial epoch, and (C) the unit directional L n in accordance with Equation (44). 
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The subroutine 132 then performs a first round of block 
back Cholesky substitution in step 154 to compute the vector 
y in the relationship shown in Equation (45). Specifically, it 
computes (A) the N vectors y iq) according to Equation (47), 
and (B) the vector Y„ according to Equation (48). 5 

In step 156, subroutine 132 performs a second round of 
block back substitution in accordance with Equation (49) to 
compute the vector x of the matrix system (40). Thus, it 
computes (A) the N vectors S0 (t?) according to Equation 
(50), and (B) the vector n according to Equation (51). The to 
flexure state value 8f and the integer ambiguities n 25(0/34/ 

35 (m> are obtained from the computed vector n in accordance 
with Equation (38). 

At step 158, subroutine 132 determines if the flexure state 
value 5f, the integer ambiguities n 25(0/34/35(m) , and the 15 
vectors 50 c<?) computed from the current iteration have 
converged to within a desired tolerance. If so, then the 
initialization process is over and the computed values n 25(0/ 

34/3 5 (m) from this last iteration represent resolution of the 
integer ambiguities are then stored in storage area 127 and 20 
used by routine 125 at step 166. As indicated earlier, routine 
125 is for computing the actual attitude determinations for 
airplane (or spacecraft) 21. 

However, if convergence has not yet been achieved, then 
subroutine 132 proceeds on to step 160. At this step, 25 
subroutine 160 computes the N skew matrices © x(<7) (i.e. the 
q th of the matrices 0* (1) through 0 x(A ° computed for the N 
measurement epochs of the initialization period) according 
to Equation (25). The subroutine 132 then computes at step 
162 the N correctional matrices 8A (<?) (i.e. the q th of the 30 
matrices 0 x(1) through 5A (<?) computed for the N measure- 
ment epochs of the initialization period) using Equation 
(23). 

At step 164, the subroutine 125 then computes and 35 
squares off (i.e., re-orthonormalizes) the N matrices A 0 C<?) 
using Equation (22). The newly computed matrices A 0 (<7) 
from this last iteration are then used at step 148 for the next 
iteration. The subroutine 164 then continues on in the way 
just described until it determines in step 158 that the flexure ^ 
state value 5f, the integer ambiguities n 25(0/ 3 4/ 3 5(m) , and the 
vectors 50 (<?) computed from the current iteration have 
converged. The routine 132 then generates a status report in 
step 165 by setting an integer valid/invalid flag for each of 
the integer ambiguity resolution values n 25(l)/34/35(m) com- 45 
puted in step 156 which indicates that these computed values 
and the phase measurements 0 25(l - )/34 and $25(0/3 5(m) which 
correspond to each of these values may be used for attitude 
computations. 

Once this is done, the attitude computation routine 125 50 
computes attitude determinations for airplane (or spacecraft) 

21 at each epoch after the initialization period. The operation 
of the attitude determination routine 125 in conjunction with 
the integrity checking and integer hand off routine 126 is 
shown in FIG. 10. 55 

First, at step 166, routine 125 computes at each iteration 
the differential phase residuals 5r 24(f)/34/35(m) and the sensi- 
tivity matrix H in accordance with Equations (65) and (66). 
Thus, the differential phase residuals Sr 24(l . )/34/35(m) are com- 
puted from (A) the integer ambiguities n 25(0/34/35(m) com- 60 
puted by routines 123 or 124 which have been flagged as 
valid in step 165 of routine 132 or in step 194 of routine 126 
(as described shortly, this means that these values can be 
used for attitude determination), (B) the measurements 
$25(0/34 and $25(0/35 (m) mad e b y rout i ne 121 at each epoch 65 
which are associated with the valid integer ambiguities 
n 25(o/34/35(m> and m tbem selves flagged as valid in step 190 
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(as described shortly, this means that the GPS signals 25(0 
associated with these measurements have sufficient signal 
strength), (C) the unit directional vectors_s 24(0/34 computed 
by routine 122, (D) the position vectors k 35(m) , and (E) the 
attitude matrix A 0 computed by routine 125 from the pre- 
vious iteration. Furthermore, the sensitivity matrix H is 
computed from (A) the unit directional vectors s 24(0/34 
computed by routine 122, (B) the matrix K 35(m> * and (C) the 
attitude matrix Aq computed by routine 125 from the pre- 
vious epoch. 

Then, at step 168, the subroutine 125 computes the vector 
50 and the flexure state correction value 88f by minimizing 
the linearized cost function represented by the matrix system 
(65). And, at step 170, the routine 125 computes the matrix 
0* using Equation (25) and the flexure state value 5f using 
Equation (60). The routine 125 then computes the correc- 
tional matrix 5A using Equation (23) at step 172. At step 
174, the routine 125 then computes the matrix A using 
Equation (22). 

In step 176, subroutine 125 determines if the flexure state 
value 8f and the attitude matrix A computed from the current 
iteration have converged to within a desired tolerance. If so, 
then the process proceeds to step 178 of the integrity 
checking and integer hand off routine 126. 

However, if convergence has not yet been achieved, then 
routine 125 returns to step 166. At this step, subroutine 125 
computes the differential phase residuals 5r 24(0/34/35(m) and 
the sensitivity matrix H in accordance with Equations (64) 
and (66) using the just computed attitude matrix A from the 
last iteration as the estimate A 0 for the next iteration. The 
routine 125 then iteratively computes the flexure state value 
8f and the attitude matrix A in the way just described until 
it determines in step 174 that the flexure state value Sf and 
the attitude matrix A from the current iteration have con- 
verged to within a desired tolerance. 

Once the flexure state value 5f and the attitude matrix A 
have converged, the integrity checking and integer handoff 
routine 126 determines if an integer slip has occurred. This 
is done as follows. 

At step 178, the routine 126 computes the RMS value of 
the differential phase residuals 5r 24(l)/34/35(m) computed in 
step 166. In step 180, the routine 126 determines if the RMS 
value calculated in step 178 exceeds a predefined threshold 
Trms- If the RMS value does not exceed the threshold T^^ 
then an integer slip has not occurred and the routine 126 
proceeds to step 186. However, if the RMS value exceeds 
the threshold then an integer slip is suspected and 

routine 126 proceeds to step 182. 

At step 182, the routine 126 examines individually the 
absolute value of each residual 5r 24(£)/34/35(m) and determines 
which of these has the highest absolute value. For those 
residuals 5r 24(0/34/35(m) that do exceed the threshold Tres, 
the routine 126 updates in step 184 the status report by 
re-setting the integer valid/invalid flag for the integer ambi- 
guity resolution value n 24(0/34/35(m) associated with the 
residual Sr 24(0/34/35(m) which had the highest absolute value. 
The reset integer valid/invalid flag indicates that this reso- 
lution value n 24(0/34/35(m) and the phase measurements 
$25(0/34 and $25 (/)/35 (m) which correspond to it can no 
longer be used for attitude computations. In other words, for 
the phase residual §r 24(0/34/35(m) with the highest absolute 
value, an integer slip has occurred in the integer portion $ in/ 
of one of the phase measurements $25(0/34 or $25(0/35 (m) 
associated with that residual. Thus, the phase measurements 
$25(0/34 or $25(o/35(m) associated with that residual cannot 
be used until a new integer ambiguity resolution value 
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n 24(o/34/35(m> is computed which takes into account the 
amount of the integer slip. 

In this case, the routine 125 then returns to step 166 and 
proceeds to compute the attitude matrix A for the next 
iteration. This is done without the integer ambiguity reso- 5 
lution values n 24(0/34/35(m) designated as being invalid in the 
status report and without the phase measurements $25(0/34 
and 3>25(o/35(m) associated with these resolution values. This 
process is iteratively repeated until all of the phase mea- 
surements $ 25(i)/3 4 and di 25(0/35 (m) with integer slips have 10 
been weeded out. 

However, if in step 180 the RMS value of the phase 
residuals §r 24(0/34/35(m) does not exceed the threshold T^^, 
then an integer slip has not occurred and the routine 126 
proceeds to step 186. At step 186, the routine 126 determines 15 
the signal strength of all of the GPS signals 25(1 )-(/') 
currently in view including those that just came into view. 
This is done by using the signal strength information pro- 
vided by the signal processing block 111. 

At step 188, the routine 126 determines if the signal 20 
strength for each of the GPS signals 25(1 )-(/) computed in 
step 186 exceeds a threshold value T str . If it does, then in 
step 190 the status report includes a phase valid/invalid flag 
for each of these signals which indicates that the phase 
measurements O 25(0/34 and O 25(0/35(m) for each of the 25 
signals can be used by routine 125 in step 166 for later 
computation of the attitude matrix A and can be used by 
routine 126 in step 192 for computing a new integer ambi- 
guity resolution value n 24(0/34/35(m) . However, if the com- 
puted signal strength does not exceed the threshold 
then in step 191 the phase valid/invalid flag for these signals 
indicates that the phase measurements <E> 25(i)/34 anc * $ 2 5(/)/ 

3 5(m> cannot be used by routine 125 in step 166 nor by 
routine 126 in step 192. 

35 

At step 192, the routine 126 then computes a new integer 
ambiguity resolution value n 24(0/34/35(m) if (A) the integer 
valid/invalid flag associated with that integer ambiguity 
resolution value n 24(r . )/34/35(m) indicates that the currently 
stored value is invalid (in the case of a newly acquired GPS 4Q 
signal 25(i) this flag means that the integer ambiguities 
n 24(o/34/35(m) associated with that signal have not yet been 
resolved), and (B) the phase valid/invalid flag indicates that 
the phase measurements 9? 2 g(- £ -y 34 and ^^25(i)/35(wz) needed to 
compute the new integer ambiguity resolution value n 24C0/ 45 
34/35 (m) are valid (i.e. the corresponding GPS signal 25(0 has 
enough signal strength). This is done in accordance with 
Equation (20) by rounding off the result to an integer. 

Thus, routine 126 computes each of these values using 
(A) the attitude matrix A calculated from the final iteration 50 
only with the integer ambiguity resolution values n 24(0/34/ 
35(m) designated as being valid by the integer valid/invalid 
flags and their corresponding phase measurements <E> 2 5(o/34 
and d> 25 (i)/35 (m) , and (B) the current phase measurements 
^25(0/34 and <b 25( j)/ 35 ( m) associated with the integer ambi- 55 
guity resolution values n 24(0/34/35(m) that are designated as 
invalid by the integer valid/invalid flags (These phase mea- 
surements will have an integer wavelength portion 
offset by the amount of the integer slip). 

Then, in step 194, the routine 126 sets the valid/invalid 60 
flag to indicate that the newly computed integer ambiguity 
resolution values n 25(om/35(m) in step 192 are valid. Then, 
these newly computed integer ambiguity resolution values 
n 25(o/34/35(m) together with the phase measurements $25(0/34 
and <I>2S(f)/35(m) that correspond to these values may be used 65 
by routine 125 in computing the matrix A in the way 
described earlier. 
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And finally, in step 196 the routine 126 outputs the final 
attitude matrix A computed for the current epoch as the 
complete attitude solution for that epoch. 

Referring again to FIG. 8, the synthesizer 116 and the 
reference oscillator 115 are coupled together. The synthe- 
sizer 116 uses the reference frequency signal output by the 
oscillator 115 for generating a clock signal. 

The clock signal generated by the synthesizer 116 is 
received by the signal processing stage 114 and the CPU 
118. Since the CPU 118 and the signal processing stage 114 
operate based on the same clock source, the carrier phase 
measurements made for each of the GPS signals 25(1)— (7) 
are coherent (i.e. made at the same time) with respect to each 
other. 

Many of the individual elements of the components of 
system 20 are known in the art. In fact, many are found in 
commercially available products. 

Specifically, the GPS antennas 34 and 35(l)-(3) are of the 
type commonly known as standard hemispherical microstrip 
patch antennas. 

The signal receiving stages 113(l)-(4), the signal pro- 
cessing stages 114(l)-(4), the synthesizer 115, the oscillator 
117, the signal processing routine 120, and the carrier phase 
measuring routine 121 may be of the type commonly found 
in a Trimble 4000 Series GPS receiver. These are described 
in U.S. Pat. No. 4,754,465 which is expressly incorporated 
by reference herein. 

Furthermore, although the invention is described with 
reference to a few specific embodiments, the description is 
illustrative of the invention and is not to be construed as 
limiting the invention. For example, the baseline vectors 
b 35(m) for spacecraft can also be computed as described in 
Cohen, “Attitude Determination Using GPS”, Ph.D. Thesis, 
Stanford University, December 1992 and in Cohen and 
Parkinson, “Integer Ambiguity Resolution of the GPS Car- 
rier for Spacecraft Attitude Determination”, Advances in the 
Astronautical Sciences, \A th Annual AAS Guidance and 
Control Conference, Keystone, Colo., February 1992, which 
are both hereby expressly incorporated herein. And, the 
initialization subroutines 131 and 132 as well as the routines 
125 and 126 can be utilized for spacecraft as well. 

Specifically, in the case of spacecraft, the platform flexure 
of the spacecraft can be assumed to be negligible. Thus, the 
flexure state value 5f and the vectors f (9) can be removed 
from matrix systems (29), (33), and (40) so that only the 
vectors 80 (?) and n are computed at each iteration of the 
initialization process. 

Furthermore, one skilled in the art would readily under- 
stand that the invention is applicable to attitude determina- 
tions for any land, sea, air, or space vehicle. Indeed, various 
modifications may occur to those skilled in the art without 
departing from the true spirit and scope of the invention as 
defined by the appended claims. 

What is claimed is: 

1. A GPS attitude system for determining the attitude of 
a moving vehicle from a plurality of GPS carrier signals, the 
GPS attitude system comprising: 

four antennas mounted to the moving vehicle, each of the 
antennas receiving the GPS carrier signals; 

signal receiving means for receiving the GPS carrier 
signals from each of the antennas; 

phase measuring means for making phase measurements 
for the GPS carrier signals received from each of the 
antennas at a plurality of epochs during an initialization 
period and at a plurality of epochs after the initializa- 
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tion period, wherein the vehicle is moving while the 
phase measurements are made during and after the 
initialization period, and wherein there is an integer 
wavelength ambiguity associated with the phase mea- 
surements made for each of the GPS carrier signals 5 
received from a first one and a second one of the 
antennas, there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
of the GPS carrier signals received from the first one 
and a third one of the antennas, and there is an integer 1Q 
wavelength ambiguity associated with the phase mea- 
surements made for each of the GPS carrier signals 
received from the first one and a fourth one of the 
antennas; 

baseline computation means for computing baseline posi- 
tions for the second, third, and fourth ones of the 1 
antennas with respect to the first one of the antennas at 
an initial one of the epochs during the initialization 
period and the changes in the baseline positions at 
subsequent ones of the epochs during the initialization 
period in response to the phase measurements made at 20 
the epochs during the initialization period and based on 
the constraint relationship at each of the subsequent 
ones of the epochs during the initialization period 

Ab z T b y +Ab y T b=-Ab y T Ab z , 25 

where b y is the baseline position of the y th of the second, 
third, and fourth ones of the antennas, b z is the baseline 
position of the z lh of the second, third, and fourth ones 
of the antennas, Ab y is the change in the baseline 30 
position of the y th of the second, third, and fourth ones 
of the antennas at a respective one of the subsequent 
epochs during the initialization period, Ab z is the 
change in the baseline position of the z th of the second, 
third, and fourth ones of the antennas at the respective 35 
one of the subsequent epochs during the initialization 
period, and Ab y r Ab z is non-zero; 
initial guess computation means for computing an initial 
guess of the attitude of the moving vehicle at each of 
the epochs during the initialization period in response 40 
to the computed baseline positions and the computed 
changes in the baseline positions; 
ambiguity resolution means for resolving the integer 
wavelength ambiguities in response to the computed 
initial guesses and the phase measurements made at the 45 
epochs during the initialization period; and 
attitude computation means for computing the attitude of 
the moving vehicle at each respective one of the epochs 
after the initialization period in response to the resolved 
integer wavelength ambiguities and the phase measure- 50 
ments made at the respective one of the epochs after the 
initialization period. 

2. The GPS system of claim 1 wherein: 
the ambiguity resolution means resolves the integer wave- 55 
length ambiguities in accordance with a predefined set 
of simultaneous equations that relate (A) the integer 
wavelength ambiguities, (B) the phase measurements 
made at the epochs during the initialization period, (C) 
estimates of the attitudes of the moving vehicle at the 
epochs during the initialization period, and (D) correc- 
tions to the estimates; 

the number of GPS carrier signals and epochs during the 
initialization period is such that the predefined set of 
simultaneous equations is overdetermined; 65 

the ambiguity resolution means iteratively resolves the 
integer wavelength ambiguities by (A) computing the 
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integer ambiguities and the corrections in a plurality of 
iterations in accordance with the predefined set of 
equations, (B) in an initial one of the iterations, using 
the initial guesses as the estimates, and (C) in each 
subsequent one of the iterations, using as the estimates 
the estimates used in a previous one of the iterations 
corrected with the corrections computed in the previous 
one of the iterations. 

3. The GPS system of claim 2 wherein the ambiguity 
resolution means resolves the integer wavelength ambigu- 
ities in accordance with block cholesky factorization of the 
predefined set of equations. 

4. The GPS system of claim 1 wherein the ambiguity 
resolution means resolves the integer wavelength ambigu- 
ities by simultaneously resolving the integer wavelength 
ambiguities and determining the state of flexure of the 
moving vehicle during the initialization period. 

5. The GPS system of claim 4 wherein: 

the ambiguity resolution means simultaneously resolves 
the integer ambiguities and determines the state of 
flexure of the moving vehicle in accordance with a 
predefined set of simultaneous equations that relate (A) 
the integer wavelength ambiguities, (B) the phase mea- 
surements made at the epochs during the initialization 
period, (C) estimates of the attitudes of the moving 
vehicle at the epochs during the initialization period, 
(D) corrections to the estimates, and (E) the state of 
flexure; 

the number of GPS carrier signals and epochs during the 
initialization period is such that the predefined set of 
simultaneous equations is overdetermined; 

the ambiguity resolution means iteratively and simulta- 
neously resolves the integer wavelength ambiguities 
and determines the state of flexure of the moving 
vehicle by (A) computing the integer ambiguities, the 
corrections, and the state of flexure in a plurality of 
iterations in accordance with the predefined set of 
equations, (B) in an initial one of the iterations, using 
the initial guesses as the estimates, and (C) in each 
subsequent one of the iterations, using as the estimates 
the estimates used in a previous one of the iterations 
corrected by the corrections computed in the previous 
one of the iterations. 

6. The GPS system of claim 5 wherein the ambiguity 
resolution means resolves the integer wavelength ambigu- 
ities in accordance with block cholesky factorization of the 
predefined set of equations. 

7. A method of determining the attitude of a moving 
vehicle from a plurality of GPS carrier signals, the method 
comprising the steps of: 

receiving the GPS carrier signals with each of four 
antennas mounted to the moving vehicle; 

receiving the GPS carrier signals from each of the anten- 
nas; 

making phase measurements for the GPS carrier signals 
received from each of the antennas at a plurality of 
epochs during an initialization period and at a plurality 
of epochs after the initialization period, wherein the 
vehicle is moving while the phase measurements are 
made during and after the initialization period, and 
wherein there is an integer wavelength ambiguity asso- 
ciated with the phase measurements made for each of 
the GPS carrier signals received from a first one and a 
second one of the antennas, there is an integer wave- 
length ambiguity associated with the phase measure- 
ments made for each of the GPS carrier signals 
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received from the first one and a third one of the 
antennas, and there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
of the GPS carrier signals received from the first one 
and a fourth one of the antennas; 5 

computing baseline positions for the second, third, and 
fourth ones of the antennas with respect to the first one 
of the antennas at an initial one of the epochs during the 
initialization period and the changes in the baseline 
positions at subsequent ones of the epochs during the 10 
initialization period in response to the phase measure- 
ments made at the epochs during the initialization 
period and based on the constraint relationship at each 
of the subsequent ones of the epochs during the ini- 
tialization period 

Ab z T b y +Ab y T b z =-Ab y T Ab z , 

where b y is the baseline position of the y !h of the second, 
third, and fourth ones of the antennas, b z is the baseline 
position of the z th of the second, third, and fourth ones 20 
of the antennas, Ab^ is the change in the baseline 
position of the y th of the second, third, and fourth ones 
of the antennas at a respective one of the subsequent 
epochs during the initialization period, Ab z is the 
change in the baseline position of the z th of the second, 25 
third, and fourth ones of the antennas at the respective 
one of the subsequent epochs during the initialization 
period, and Ab y r Ab z is non-zero; 
computing an initial guess of the attitude of the moving 
vehicle at each of the epochs during the initialization 30 
period in response to the computed baseline positions 
and the computed changes in the baseline positions; 
resolving the integer wavelength ambiguities in response 
to the computed initial guesses and the phase measure- 
ments made at the epochs during the initialization 35 
period; and 

computing the attitude of the moving vehicle at each 
respective one of the epochs after the initialization 
period in response to the resolved integer wavelength 
ambiguities and the phase measurements made at the 40 
respective one of the epochs after the initialization 
period. 

8. The method of claim 7 wherein: 

the integer wavelength ambiguities are resolved in the 45 
ambiguity resolving step in accordance with a pre- 
defined set of simultaneous equations that relate (A) the 
integer wavelength ambiguities, (B) the phase measure- 
ments made at the epochs during the initialization 
period, (C) estimates of the attitudes of the moving 5Q 
vehicle at the epochs during the initialization period, 
and (D) corrections to the estimates; 
the number of GPS carrier signals and epochs during the 
initialization period is such that the predefined set of 
simultaneous equations is overdetermined; 55 

the integer wavelength ambiguities are iteratively 
resolved in accordance with the predefined set of 
equations in the ambiguity resolving step by (A) com- 
puting the integer ambiguities and the corrections in a 
plurality of iterations, (B) in an initial one of the 60 
iterations, using the initial guesses as the estimates, and 
(C) in each subsequent one of the iterations, using as 
the estimates the estimates used in a previous one of the 
iterations corrected with the corrections computed in 
the previous one of the iterations. 65 

9. The method of claim 8 wherein the integer wavelength 
ambiguities are iteratively resolved in the ambiguity resolv- 
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ing step in accordance with block cholesky factorization of 
the predefined set of equations. 

10. The method of claim 7 wherein the integer wavelength 
ambiguities are resolved in the ambiguity resolving step by 
simultaneously resolving the integer wavelength ambigu- 
ities and determining the state of flexure of the moving 
vehicle during the initialization period. 

11. The method of claim 10 wherein: 

the integer wavelength ambiguities are resolved and the 
state of flexure of the moving vehicle is determined in 
the ambiguity resolving step in accordance with a 
predefined set of simultaneous equations that relate (A) 
the integer wavelength ambiguities, (B) the phase mea- 
surements made at the epochs during the initialization 
period, (C) estimates of the attitudes of the moving 
vehicle at the epochs during the initialization period, 
(D) corrections to the estimates, and (E) the state of 
flexure; 

the number of GPS carrier signals and epochs during the 
initialization period is such that the predefined set of 
simultaneous equations is overdetermined; 

the integer wavelength ambiguities are resolved and the 
state of flexure of the moving vehicle is determined 
iteratively and simultaneously in the ambiguity resolv- 
ing step by (A) computing the integer ambiguities, the 
corrections, and the state of flexure in a plurality of 
iterations in accordance with the predefined set of 
equations, (B) in an initial one of the iterations, using 
the initial guesses as the estimates, and (C) in each 
subsequent one of the iterations, using as the estimates 
the estimates used in a previous one of the iterations 
corrected by the corrections computed in the previous 
one of the iterations. 

12. The method of claim 11 wherein the integer wave- 
length ambiguities are resolved and the state of flexure of the 
moving vehicle is determined iteratively and simultaneously 
in the ambiguity resolving step in accordance with block 
cholesky factorization of the predefined set of equations. 

13. A GPS attitude system for determining the attitude of 
a moving vehicle from a plurality of GPS carrier signals, the 
GPS attitude system comprising: 

four antennas mounted to the moving vehicle, each of the 
antennas receiving the GPS carrier signals; 

signal receiving means for receiving the GPS carrier 
signals from each of the antennas; 

phase measuring means for making phase measurements 
for the GPS carrier signals received from each of the 
antennas at a plurality of epochs during an initialization 
period and at a plurality of epochs after the initializa- 
tion period, wherein there is an integer wavelength 
ambiguity associated with the phase measurements 
made for each of the GPS carrier signals received from 
a first one and a second one of the antennas, there is an 
integer wavelength ambiguity associated with the 
phase measurements made for each of the GPS carrier 
signals received from the first one and a third one of the 
antennas, and there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
of the GPS carrier signals received from the first one 
and a fourth one of the antennas; 

ambiguity resolution means for simultaneously resolving 
the integer wavelength ambiguities and determining the 
state of flexure of the moving vehicle during the 
initialization period in response to the phase measure- 
ments made at the epochs during the initialization 
period; 
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attitude computation means for computing the attitude of 
the moving vehicle at each respective one of the epochs 
after the initialization period in response to the resolved 
integer wavelength ambiguities and the phase measure- 
ments made at the respective one of the epochs after the 5 
initialization period. 

14. The GPS system of claim 13 further comprising: 
initial guess computation means for computing an initial 

guess of the attitude of the moving vehicle at each of 
the epochs during the initialization period; 10 

the ambiguity resolution means simultaneously resolving 
the integer wavelength ambiguities and determining the 
state of flexure of the moving vehicle in accordance 
with a predefined set of simultaneous equations that 
relate (A) the integer wavelength ambiguities, (B) the 15 
phase measurements made at the epochs during the 
initialization period, (C) estimates of the attitudes of 
the moving vehicle at the epochs during the initializa- 
tion period, (D) corrections to the estimates, and (E) the 
state of flexure; 

20 

the number of GPS carrier signals and epochs during the 
initialization period being such that the predefined set 
of simultaneous equations is overdetermined; 
the ambiguity resolution means iteratively and simulta- 
neously resolving the integer wavelength ambiguities 25 
and determining the state of flexure of the moving 
vehicle by (A) computing the integer wavelength ambi- 
guities, the corrections, and the state of flexure in a 
plurality of iterations in accordance with the predefined 
set of simultaneous equations, (B) in an initial one of 
the iterations, using the initial guesses as the estimates, 30 
and (C) in each subsequent one of the iterations, using 
as the estimates the estimates used in a previous one of 
the iterations corrected by the corrections computed in 
the previous one of the iterations. 

15. The GPS system of claim 14 wherein the ambiguity 35 
resolution means resolves the integer wavelength ambigu- 
ities in accordance with block cholesky factorization of the 
predefined set of equations. 

16. The GPS system of claim 13 wherein the attitude 
computation means simultaneously computes the attitude 40 
and determines the state of flexure of the moving vehicle at 
each one of the epochs after the initialization period. 

17. A method of determining the attitude of a moving 

vehicle from a plurality of GPS carrier signals, the method 
comprising: 45 

receiving the GPS carrier signals with each of four 
antennas mounted to the moving vehicle; 
receiving the GPS carrier signals from each of the anten- 

n3S; 50 

making phase measurements for the GPS carrier signals 
received from each of the antennas at a plurality of 
epochs during an initialization period and at a plurality 
of epochs after the initialization period, wherein there 
is an integer wavelength ambiguity associated with the 55 
phase measurements made for each of the GPS carrier 
signals received from a first one and a second one of the 
antennas, there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
of the GPS carrier signals received from the first one 50 
and a third one of the antennas, and there is an integer 
wavelength ambiguity associated with the phase mea- 
surements made for each of the GPS carrier signals 
received from the first one and a fourth one of the 
antennas; 65 

simultaneously resolving the integer wavelength ambigu- 
ities and determining the state of flexure of the moving 
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vehicle during the initialization period in response to 
the phase measurements made at the epochs during the 
initialization period; 

computing the attitude of the moving vehicle at each 
respective one of the epochs after the initialization 
period in response to the resolved integer wavelength 
ambiguities and the phase measurements made at the 
respective one of the epochs after the initialization 
period. 

18. The GPS system of claim 17 further comprising the 
step of: 

computing an initial guess of the attitude of the moving 
vehicle at each of the epochs during the initialization 
period; 

the integer wavelength ambiguities are resolved and the 
state of flexure of the moving vehicle are determined 
simultaneously in the ambiguity resolving step in 
accordance with a predefined set of simultaneous equa- 
tions that relate (A) the integer wavelength ambiguities, 
(B) the phase measurements made at the epochs during 
the initialization period, (C) estimates of the attitudes 
of the moving vehicle at the epochs during the initial- 
ization period, (D) corrections to the estimates, and (E) 
the state of flexure; 

the number of GPS carrier signals and epochs during the 
initialization period being such that the predefined set 
of simultaneous equations is overdetermined; 

the integer wavelength ambiguities being resolved and the 
state of flexure of the moving vehicle being determined 
iteratively and simultaneously by (A) computing the 
integer wavelength ambiguities, the corrections, and 
the state of flexure in a plurality of iterations in accor- 
dance with the predefined set of simultaneous equa- 
tions, (B) in an initial one of the iterations, using the 
initial guesses as the estimates, and (C) in each subse- 
quent one of the iterations, using as the estimates the 
estimates used in a previous one of the iterations 
corrected by the corrections computed in the previous 
one of the iterations. 

19. The method of claim 18 wherein the integer wave- 
length ambiguities are resolved and the state of flexure of the 
moving vehicle is determined iteratively and simultaneously 
in the ambiguity resolving step in accordance with block 
cholesky factorization of the predefined set of equations. 

20. The method of claim 17 wherein the attitude of the 
moving vehicle is computed and the state of flexure of the 
moving vehicle is determined simultaneously at each one of 
the epochs after the initialization period in the attitude 
computing step. 

21. A GPS attitude system for determining the attitude of 
a moving vehicle from a plurality of GPS carrier signals, the 
GPS atdtude system comprising: 

a plurality of antennas mounted to the moving vehicle, 
each of the antennas receiving the GPS carrier signals; 

signal receiving means for receiving the GPS signals from 
each of the antennas; 

phase measuring means for making phase measurements 
for the GPS carrier signals received from each of the 
antennas at a plurality of epochs; 

attitude computation means for simultaneously comput- 
ing the attitude and determining the state of flexure of 
the moving vehicle at each respective one of the epochs 
in response to the phase measurements made at the 
respective one of the epochs. 

22. The GPS attitude system of claim 21 wherein: 

the plurality of antennas comprises four antennas; 
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there is an integer wavelength ambiguity associated with 
the phase measurements made for each of the GPS 
carrier signals received from a first one and a second 
one of the antennas, there is an integer wavelength 
ambiguity associated with the phase measurements 5 
made for each of the GPS carrier signals received from 
the first one and a third one of the antennas, and there 
is an integer wavelength ambiguity associated with the 
phase measurements made for each of the GPS carrier 
signals received from the first one and a fourth one of 10 
the antennas; 

the GPS attitude system further comprises ambiguity 
resolution means for resolving the integer wavelength 
ambiguities; 

the attitude computation means simultaneously computes 
the attitude and determines the state of flexure of the 
moving vehicle at each respective one of the epochs 
also in response to the resolved integer wavelength 
ambiguities. 

23. The GPS system of claim 22 wherein: 
the attitude computation means simultaneously computes 

the attitude and determines the state of flexure of the 
moving vehicle at each respective one of the epochs 
after the initialization period in accordance with a 
predefined set of simultaneous equations that relate (A) 
the resolved integer wavelength ambiguities, (B) the 
phase measurements made at the respective one of the 
epochs, (C) an estimate of the attitude of the moving 
vehicle at the respective one of the epochs, (D) a 3Q 
correction to the estimate of the attitude, (E) an esti- 
mate of the state of flexure of the moving vehicle at the 
respective one of the epochs, and (F) a correction to the 
estimate of the state of flexure; 
the attitude computation means iteratively and simulta- 35 
neously computes the attitude and determines the state 
of flexure of the moving vehicle at each respective one 
of the epochs by (A) computing the corrections to the 
estimates of the attitude and the state of flexure in a 
plurality of iterations in accordance with the predefined 40 
set of simultaneous equations, (B) in each respective 
one of the iterations, using as the estimates of the 
attitude and the state of flexure the estimates of the 
attitude and the state of flexure used in a directly 
preceding one of the iterations corrected by the correc- 45 
tions to the estimates of the attitude and state of flexure 
computed in the directly preceding one of the itera- 
tions. 

24. A method of determining the attitude of a moving 
vehicle from a plurality of GPS carrier signals, the method 50 
comprising: 

receiving the GPS carrier signals with each of a plurality 
of antennas mounted to the moving vehicle; 
receiving the GPS signals from each of the antennas; 
making phase measurements for the GPS carrier signals 
received from each of the antennas at a plurality of 
epochs; 

simultaneously computing the attitude and determining 
the state of flexure of the moving vehicle at each 
respective one of the epochs in response to the phase 
measurements made at the respective one of the 
epochs. 

25. The method of claim 24 wherein: 
the plurality of antennas comprises four antennas; 
there is an integer wavelength ambiguity associated with 

the phase measurements made for each of the GPS 
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carrier signals received from a first one and a second 
one of the antennas, there is an integer wavelength 
ambiguity associated with the phase measurements 
made for each of the GPS carrier signals received from 
the first one and a third one of the antennas, and there 
is an integer wavelength ambiguity associated with the 
phase measurements made for each of the GPS carrier 
signals received from the first one and a fourth one of 
the antennas; 

the method further comprises the step of resolving the 
integer wavelength ambiguities; 
the attitude of the moving vehicle is computed and the 
state of flexure of the moving vehicle is determined 
simultaneously at each respective one of the epochs in 
15 the attitude computing step also in response to the 
resolved integer wavelength ambiguities. 

26. The method of claim 25 wherein: 
the attitude of the moving vehicle is computed and the 

state of flexure of the moving vehicle is determined 
simultaneously at each respective one of the epochs 
after the initialization period in the attitude computing 
step in accordance with a predefined set of simulta- 
neous equations that relate (A) the resolved integer 
wavelength ambiguities, (B) the phase measurements 
made at the respective one of the epochs, (C) an 
estimate of the attitude of the moving vehicle at the 
respective one of the epochs, (D) a correction to the 
estimate of the attitude, (E) an estimate of the state of 
flexure of the moving vehicle at the respective one of 
the epochs, and (F) a correction to the estimate of the 
state of flexure; 

the attitude of the moving vehicle is computed and the 
state of flexure of the moving vehicle is determined 
iteratively and simultaneously at each respective one of 
the epochs after the initialization period in the attitude 
computing step by (A) computing the corrections to the 
estimates of the attitude and the state of flexure in a 
plurality of iterations in accordance with the predefined 
set of simultaneous equations, (B) in each respective 
one of the iterations, using as the estimates of the 
attitude and the state of flexure the estimates of the 
attitude and the state of flexure used in a directly 
preceding one of the iterations corrected by the correc- 
tions to the estimates of the attitude and state of flexure 
computed in the directly preceding one of the itera- 
tions. 

27. A GPS attitude system for determining the attitude of 
a moving vehicle from a plurality of GPS carrier signals, the 
GPS attitude system comprising: 

four antennas mounted to the moving vehicle, each of the 
antennas receiving the GPS carrier signals; 
signal receiving means for receiving the GPS signals from 
each of the antennas; 

55 phase measuring means for making phase measurements 
for the GPS carrier signals received from each of the 
antennas at a plurality of epochs during an initialization 
period and at a plurality of epochs after the initializa- 
tion period, wherein there is an integer wavelength 
ambiguity associated with the phase measurements 
made for each of the GPS carrier signals received from 
a first one and a second one of the antennas, there is an 
integer wavelength ambiguity associated with the 
phase measurements made for each of the GPS carrier 
65 signals received from the first one and a third one of the 
antennas, and there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
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of the GPS carrier signals received from the first one 
and a fourth one of the antennas; 
ambiguity resolution means for (A) resolving the integer 
wavelength ambiguities, and (B) generating a status 
report indicating that the resolved integer ambiguities 5 
are valid; 

integrity checking means for (A) determining when an 
erroneous integer wavelength change has occurred in 
the phase measurement made for a particular one of the 
GPS carrier signals received by a particular one of the 10 
antennas at a particular one of the epochs after the 
initialization period, the erroneous integer wavelength 
change creating one or more new integer wavelength 
ambiguities associated with the phase measurements 
made for the particular one of the GPS carrier signals 15 
received by the particular one of the antennas at suc- 
ceeding ones of the epochs after the initialization 
period, and (B) updating the status report to indicate 
that each one of the resolved integer wavelength ambi- 
guities associated with the phase measurement made 20 
for the particular one of the GPS carrier signals 
received by the particular one of the antennas at the 
particular one of the epochs is no longer valid; 
attitude computation means for computing the attitude of 25 
the moving vehicle at each respective one of the epochs 
after the initialization period in response to the resolved 
integer wavelength ambiguities designated as being 
valid in the status report and the phase measurements 
made at the respective one of the epochs after the 3Q 
initialization period associated with the resolved inte- 
ger ambiguities designated as being valid in the status 
report; 

the integrity checking means also (C) resolving the new 
integer wavelength ambiguities in response to the atti- 35 
tude of the moving vehicle computed at an initial one 
of the succeeding ones of the epochs after the initial- 
ization period and the phase measurements made at the 
initial one of the succeeding ones of the epochs after the 
initialization period associated with the new integer 40 
wavelength ambiguities, and (D) updating the status 
report to indicate that the resolved new integer wave- 
length ambiguities are valid. 

28. A method of determining the attitude of a moving 
vehicle from a plurality of GPS carrier signals, the method 45 
comprising: 

receiving the GPS carrier signals with each of four 
antennas mounted to the moving vehicle; 
receiving the GPS signals from each of the antennas; 
making phase measurements for the GPS carrier signals 50 
received from each of the antennas at a plurality of 
epochs during an initialization period and at a plurality 
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of epochs after the initialization period, wherein there 
is an integer wavelength ambiguity associated with the 
phase measurements made for each of the GPS carrier 
signals received from a first one and a second one of the 
antennas, there is an integer wavelength ambiguity 
associated with the phase measurements made for each 
of the GPS carrier signals received from the first one 
and a third one of the antennas, and there is an integer 
wavelength ambiguity associated with the phase mea- 
surements made for each of the GPS carrier signals 
received from the first one and a fourth one of the 
antennas; 

resolving the integer wavelength ambiguities; 

generating a status report indicating that the resolved 
integer ambiguities are valid; 

determining when an erroneous integer wavelength 
change has occurred in the phase measurement made 
for a particular one of the GPS carrier signals received 
by a particular one of the antennas at a particular one 
of the epochs after the initialization period, the erro- 
neous integer wavelength change creating one or more 
new integer wavelength ambiguities associated with 
the phase measurements made for the particular one of 
the GPS carrier signals received by the particular one of 
the antennas at succeeding ones of the epochs after the 
initialization period; 

updating the status report to indicate that each one of the 
resolved integer wavelength ambiguities associated 
with the phase measurement made for the particular 
one of the GPS carrier signals received by the particular 
one of the antennas at the particular one of the epochs 
is no longer valid; 

computing the attitude of the moving vehicle at each 
respective one of the epochs after the initialization 
period in response to the resolved integer wavelength 
ambiguities designated as being valid in the status 
report and the phase measurements made at the respec- 
tive one of the epochs after the initialization period 
associated with the resolved integer ambiguities; 

resolving the new integer wavelength ambiguities in 
response to the attitude of the moving vehicle com- 
puted at an initial one of the succeeding ones of the 
epochs after the initialization period and the phase 
measurements made at the initial one of the succeeding 
ones of the epochs after the initialization period asso- 
ciated with the new integer wavelength ambiguities; 
and 

updating the status report to indicate that the resolved new 
integer wavelength ambiguities are valid. 



